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cosmic string loops, which simulations suggest would be far more numerous than 
long (horizon-sized) strings. Using simple models of the loop population and minimal 
assumptions about the lensing cross section per loop, we estimate the optical depth 
to lensing and show that extant radio surveys such as CLASS have already ruled out 
a portion of the cosmic string model parameter space. Future radio interferometers, 
such as LOFAR and especially SKA, may constrain Gfi/c^ < 10~^ in some regions 
of parameter space, outperforming current constraints from pulsar timing and the 
CMB by up to two orders of magnitude. This method relies on direct detections 
of cosmic strings, and so is less sensitive to the theoretical uncertainties in string 
network evolution that weaken other constraints. 
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I. INTRODUCTION 



Cosmic strings are linear, relativistic objects whose existence is predicted by a number 
of extensions to the Standard Model of particle physics When first proposed, 

they were intensely studied as a possible mechanism for cosmological structure formation. 
This scenario is now known to be inconsistent with a variety of cosmological constraints 
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I2I . [13I . Il^ , but recently cosmic strings have attracted renewed interest 



thanks to the postulation of new types and new production mechanisms. The current focus 
has shifted from their role in structure formation to the prospects for detecting them if they 



are present 
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Here we describe a method to detect cosmic string loops using strong gravitational lensing 
of compact radio sources. We propose searching for loops because simulations suggest there 
could be > 10^ per horizon volume in contrast to only ~ 40 long (horizon-sized) strings 
Compact radio sources (hereafter CRSs) are an ideal source population since their point- 
like nature on angular scales < 100 mas makes lensed images easy to identify. Furthermore, 
numerous CRSs have been observed by radio surveys such as the Cosmic Lens All-Sky Survey 
(CLASS) and are expected to be observed by proposed surveys such as the Low Frequency 
Array (LOFAR) and the Square Kilometer Array (SKA). In this paper, we estimate the 
expected number of loop lensing events for these three surveys as a function of cosmic string 
model parameters, and find that for significant regions of parameter space our method can 
outperform existing constraints. 

An important feature of our method is that it would constitute a direct detection rather 
than an inference based on assumed or simulated properties of a cosmic string network. 
By contrast, the pulsar timing constraint (the most stringent constraint to date) requires 
knowledge of the rate of emission of gravitational waves by oscillating strings and of prop- 
erties of the loop population. As we will discuss below, the signatures of lensing by a string 
loop are distinctively different from those of lensing by a galaxy (which will generally have 
a similar image separation in the arcsecond regime). Follow-up observations should be able 
to quickly determine if putative lensing events are due to cosmic strings or to more conven- 
tional astronomical objects. In past gravitational lensing surveys, all lensing events have 
been associated with conventional lensing objects such as galaxies or clusters after follow- 
up observations. We rely on the reasonable expectation that sufficient follow-up is always 
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performed; a constraint will be obtained if all lensing events are accounted for by standard 
lenses and are therefore not cosmic string lensing candidates. If after deep imaging a conven- 
tional lensing object was not observed, the lens would be classified as "dark" and modelling 
of the gravitational potential would be required to rule it out as a cosmic string. As of yet, 
no confirmed dark lenses have been found. 

Either a detection of, or a constraint on, cosmic strings would be important for under- 
standing physics beyond the Standard Model. Cosmic strings can be accommodated in 
many superstring models, and indeed are predicted to be present in brane inflation scenar- 
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25l |. The wide range of string tensions and the novel features of 



cosmic strings predicted by these models present a challenge to observation. By tightening 
constraints on cosmic strings, we are indirectly testing these models and learning about 
physics at high energies and early times in the history of the universe. 

This paper is organized as follows. In Section [TTl we describe our assumptions about the 
properties of the loop lenses, discuss our population models, and describe the telescopes and 
surveys proposed for the observation. We present the results of the calculation in Section 
llllt giving the expected number of lensing events as a function of parameters of the string 
model and comparing the constraint obtained to those presently available. We conclude in 
Section IIVI and discuss areas in which future work is needed. 



II. METHODOLOGY 



We now describe our calculation of the expected number of observable string loop lensing 
events. We do this by combining estimates of the area of sky that can be strongly lensed 
by loops, the resolution and dynamic range of the instrument employed, and the source 
population abundance, which is a function of instrument and survey parameters as well as 
of the intrinsic abundance. 

To estimate of the fraction of the sky that is strongly lensed by the loop population we 
must confront the significant theoretical uncertainties in both (i) the lensing characteristics 
of individual loops, and (ii) the distribution of lengths in the loop population as a function of 
redshift. To deal with (i) we take a statistical approach: we assume that each loop can form 
multiple images of all sources in a patch of sky of angular area 116^, with image separations 
of at least 6e, where 6'^; is the Einstein radius of a Schwarzschild lens with the same mass 
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as the loop. 

Realistic loops will be of highly irregular shape and can produce complicated lensing 



patterns, which have been considered in detail elsewhere [26|, |27|, |28|, |29|]. However, we 
expect that the lensing cross section, when averaged over a population of loops of fixed 
mass, will be close to the characteristic area set by the Schwarzschild lens. We justify our 
approximation in Section Hi Al by comparing it to the exact results in cases where loop lensing 
can be studied analytically. To model the loop population as required by (ii), we use two 
simple models, which we discuss in Section IIIBt with the string tension fi a free parameter 
in both. The first is a powerlaw distribution of loop lengths, defined by a spectral index 
7, the loop density parameter Qioop, and a cutoff length L^. The second model is derived 
from the one-scale model, a semi-analytic model of string networks widely discussed in the 
literature. This model has a number of parameters, but in addition to fi we consider only a 
(which sets the size of newly created loops) to be free. 

We then turn to the source population and instruments. In Section III CI we argue that 
compact radio sources make ideal source candidates for investigating lensing by loops, and we 
discuss some of their features. In Section Hi D| we describe both extant surveys (CLASS) and 
proposed instruments (LOFAR and SKA) that will observe these objects. Our discussion 
focuses on the parameters of these surveys that are most relevant to a prediction of the 
number of string loop lensing events expected in each survey, namely the angular resolution 
and sensitivity. The angular resolution sets the minimum separation of images which may 
be detected, and hence the minimum string tension of observable loops. The sensitivity of 
the instrument influences the number of potential targets for lensing, since a more sensitive 
instrument will pick up fainter or more distant sources. 

A. Loops as lenses 

We rely on the approximation, justified below, that a loop of mass Mioop produces multiple 
images of background sources situated in a patch of sky of area ttO"^, and that the images 
are separated by at least 6^, where 6e is the Einstein radius of a Schwarzschild lens having 
the same mass as the loop. The Einstein radius for a loop of mass Mioop is given by 



= \l D,D 
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where Dl, Ds and Dls are the angular diameter distances between the observer and lens, 
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observer and source, and lens and source, respectively 

To justify this approximation we compare it to the sky area lensed by a circular loop 
of radius R and angular radius Oi^op = R/D^ situated in a plane perpendicular to the line 
of sight. Light rays passing outside the loop are deflected as if by a Schwarzschild lens of 
the same mass as the loop, while those passing through the loop are not deflected. A pure 
Schwarzschild lens always produces two images of a source at angles 6± from the lens, given 

where 6s is the (unlensed) angular position of the source. For a loop, any images at ^ < Oioop 
are spurious, since the corresponding light rays must have passed through the loop, and will 
not be observed. The unlensed source itself will be visible through the loop if 9s < Oioop- 

The region of {6s, dioop) parameter space in which multiple images form is shown in 
Figure [H While multiple images are always created by a Schwarzschild lens, the "+" and 
"— " images are produced with different fluxes, and if the flux ratio is too large then the pair 
will not be observed as a lens. Imposing a conservative flux ratio < 10:1 for our calculations 
means that the "— " image will only be included if 6s < 1.125 6e, which cuts off the multiple 
lensing region for the "— " image. The source and "+" image have a flux ratio that is nearly 
unity, which enhances the prospects for lens identification when both are present. When 
6ioop ^ Oe there is always at least one pair of images separated by > 29 e- For 9ioop > Oe, 
the image separations gradually decrease from 9e when 9ioop ^ 9e, and go as 9e/9ioop for 
9ioop ^ 9e- The patch of sky for which multiple images form when 9ioop > 9e has area 

7r9l (2 - [9e/9ioop]^) (3) 

which is always greater than vr6'g. One can show that the sky area ([3]) asymptotes to that 
of a long straight string as 9e/9ioop — ^ 0, corresponding to i? — oo. 

Our approximation accurately reflects the lensing cross section in the special case of a 
planar circular loop. However, the area lensed by a typical loop most likely exceeds this. 
A circular loop is the least dense arrangement of the loop's mass, and furthermore rays 
passing through the loop are not lensed at all. A realistic loop will likely be "crumpled," 
concentrating the mass and better approximating a Schwarzschild lens for a larger range of 
impact parameters. In the limit that the loop can be treated as a random walk with step 
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FIG. 1: Regions of the {Os,0ioop) parameter space where multiple images are formed by a planar 
circular loop. Regions are labeled depending on whether the it Schwarzschild images are formed, 
or whether direct images of the source (S), are visible. 

size (correlation length) Ic, then a loop of radius R, when crumpled, will fit in a sphere of 
radius CyJlcR, with C a constant of order unity. This reduces the angle Oioop subtended by 



the loop by the factor C^/TjR, which will be significant for large loops. There are other 
effects that may also increase the area of the lensed patch of sky, which we will discuss in 
Section lYl 
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B. Loop population models 



We now describe our two models of the loop length distribution, defined by a number 
density per logarithmic interval in the loop length L, Nl{L,z) = dA^(z)/dlnL. Here N 
is the number of loops in a comoving volume (c/Hq)^, where c is the speed of light and 
Hq is today's Hubble parameter. The first model is the simplest scale-free distribution, a 
powerlaw in the loop length. It is not motivated by any model of string network evolution 
that we are aware of, but is useful for understanding how the constraint depends on basic 
features of the loop length spectrum and may describe other lens populations. The second 
model is more realistic. It uses the loop length spectrum derived from the one-scale model 
(OSM) of cosmic string network evolution, a semi-analytic model that fits well to full-scale 
simulations. 



1. Powerlaw loop spectrum 
For the powerlaw spectrum, assuming a constant comoving density of loops, we have 

3 / r \ 7-1 



Nl(L,z) 



loopfl ' 



Pcritfi 1 7 1 



C 



L 



(4) 



where flioopfl is the present mass parameter in loops, and pcrit,o is today's critical density. 
L^, is a loop length cutoff, so we include only L < when 7 > 0, or L > when 7 < 0. If 
we instead assume a constant number of loops per particle horizon volume we obtain 



NLmz) = n 



loop' 



Peril 



c 



H{z) 

Hn 



(5) 



To better compare with the one-scale model (discussed in the following section) we will use 
this assumption (equation [5]) in all our calculations with the powerlaw model henceforth. 
For the present work we will only consider 7 < 0, as suggested by string network simu- 



lations 
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32| . The effective lower cut-off on the powerlaw loop spectrum then comes from 



the smaller of L^, and the resolution limit. In practice, Qioop and are degenerate in the sky 
fraction calculation, so for this work we fix = 1 kpc and vary Qioop only. The upper limit 
on the loop spectrum integration is chosen so that loops longer than L„^ax = (e/(l + 
which account for a fraction e of the total number of loops, are left out. In cases where the 
instrumental resolution is high enough to resolve loops all the way down to the lower cut-off 
L*, defining the upper cut-off in this way results in the sky fraction being independent of 7. 
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2. One-scale model 



The one-scale model is a simple model of string networks, supported by numerical sim- 
ulations: a detailed exposition may be found in The only length scale in this model is 
i(t), the radius of the particle horizon at a time t. Loops that form at time t have length 
a£{t), with a a constant, and subsequently shrink by emitting gravitational radiation. The 
rate of loop formation is found by requiring that a sufficient number of loops are formed to 
keep the number of long strings per particle horizon volume constant. In principle, a can 
be determined from numerical simulations, but as yet a clear consensus on its value has not 
emerged, so we take it to be a free parameter. 

By assuming a matter-dominated universe one finds 



dN 



1 + 



23.5 



dlnL 50.3a[G/i/c2 



1 



(6) 



33.3[Gm/c2]^H(, ^ 

where ino = c/Hq, and we have used the values for other one-scale model parameters in 
[4]. This distribution has a peak at Lp^ak ~ 33£i:/,)[G/i/c^], and the spectral index 7 = 
for L > Lpeak and 7 = +2 for L < Lp^ak- It has an upper cutoff at L = ai{t) = 2aiHo 
(during the matter era); this is the length at which loops are created. Using this spectrum, 
the density parameter in loops is 
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(7) 



where y = a/{Gfi/c^). These equations are derived in detail in Appendix Rl 



C. Compact Radio Sources 

Selection of a background source population has a major impact on the efficiency and 
completeness of a gravitational lensing survey. For many background sources with extended 
emission, multiple imaging can be mistaken for the intrinsic structure of the source, or 
vice-versa. In order to easily determine whether or not lensing has occurred, the best 
strategy is to optimize the survey to observe sources that are intrinsically point-like (to 
prevent the misidentification of substructure) and to use instruments with resolution high 
enough to distinguish the images produced by a typical lensing event. With this in mind. 
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we consider a source population of flat spectrum compact radio sources (CRSs), which are 
single-component sources at ~ 100 mas resolution. These sources are numerous at high 



redshifts; typically 2 ~ 1 — 2 over a few magnitudes in flux density (see, e.g., [33|, |35| ) . 
This enables us to probe a large cosmological volume. Compact radio sources are also ideal 
for the purpose of optimizing the survey resolution, as long-baseline radio interferometry 
routinely reaches sub-arcsecond resolution. Furthermore, radio sources may be significantly 
polarized, offering a further check that lensing is indeed occurring, as well as providing 
additional constraints on lens models. 

The signatures of lensing by cosmic string loops would be quite different from those 
seen in lensing by galaxies. In every instance of lensing by a galaxy, the lens would be 
identified by its presence in follow-up optical or near-infrared observations (normal elliptical 
and spiral galaxies are fairly radio quiet). In addition, lens modeling would indicate whether 
the lens has a mass density profile compatible with a galaxy, the majority of which have 



near isothermal profiles on the scales relevant to strong lensing [34] • The polarization of 
the images would also help to differentiate between lensing by galaxies and by cosmic string 
loops. In galaxy lensing, any polarized emission of the source is affected by Faraday rotation 
as it passes through the interstellar medium of the lens galaxy, resulting in (frequency- 
dependent) differences in the polarization position angles of the images; this would not 
occur if the lens was a cosmic string loop. Using these criteria it will be straightforward to 
rule out lensing by a galaxy, and careful follow-up imaging and spectroscopic observations 
will make it possible to determine if the lens is any known astrophysical object. 

D. Radio Surveys 

Having described the fraction of sky that could be lensed by string loops, and having 
selected our source population, we now discuss some instruments and surveys that might 
observe loop lensing events. Key features that will serve as inputs to the estimated number 
of lensing events are the dynamic range, resolution, and expected source counts for each 
survey. 



CLASS/JVAS. The Cosmic Lens All-Sky Survey (CLASS) |36|, l37| and the Jodrell Bank 



VLA Astrometric Survey (JVAS) 38| searched for gravitationally lensed fiat spectrum radio 



sources. Together, they targeted 16503 sources, of which 11685 have fiat spectra (spectral 
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index flatter than 0.5 between 1.4 and 5 GHz) and flux density > 30mJy at 5 GHz. The 
statistically well deflned sample of 11 lens systems found includes only those with separations 
in excess of 0.3-arcsec and flux ratio < 10:1. 

LOFAR. The Low Frequency Radio Array (LOFAR) is expected to be operational by the 
end of 2008, initially with a maximum baseline of 100 km, and maximum resolution of ~3 
arcsec, which would only detect the largest separation events. Therefore, we concentrate on 



a proposed extension to LOFAR with a maximum baseline of 400 km [39|. Two planned 
surveys are relevant to lensing: one at 200 MHz over 250 deg^ at 0.8-arcsec resolution, and 
the other at 120 MHz over half the sky with a resolution of 1.3 arcsec 40||. The former will 
include 3x10^ sources to the 14yuJy limit and the latter 8.6x10® sources to 43/iJy. These 
source counts are expected to be dominated by steep spectrum galaxies, with ~ 5% being 
flat-spectrum 4l[ |. 



SKA. The Square Kilometre Array (SKA), planned to be fully operational by 2020, will 
consist of thousands of radio antennae with a collecting area exceeding thirty times that of 
the largest telescope ever constructed. The site and speciflcations have yet to be flnalized; 
we adopt the estimates from 42]. They discuss a possible Radio All-SKA Lens survey 
(RASKAL) of half the sky at 1.4 GHz, with 0.01-0.02-arcsec resolution. Limiting to sources 
brighter than 3/iJy there would be ~ 10^ radio sources, of which roughly 10% would be 
compact, flat-spectrum sources - the ideal target for a statistically complete lens survey. As 
with CLASS/JVAS, the median redshift is expected to be 2 > 1. 

We summarize the properties of the surveys we consider in Table [H For Extended LOFAR, 
we take the higher-resolution 240 MHz survey, and for SKA we assume the optimistic end 
of the resolution range, 0.01 arcsec. In the table, "number of sources" refers to the number 



of CRSs accessible to the survey. Based on estimates in [33| of the mean redshift of CRSs, 
for the purpose of this study we choose a redshift of 2; = 1.2 for our source population in all 
surveys. The true distribution of CRSs is still incompletely understood, especially for the 
flux limits that can be reached by future surveys. We expect that this approximation affects 
our result by less than an order of magnitude. 
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III. RESULTS 



A. Expected Lensing Events 

We estimate the number of cosmic string loop lensing events one can expect to observe 
in past and future compact radio source gravitational lensing surveys whose properties can 
be found in Table [B Our main results are summarized in Figures [2] and [31 



1. Powerlaw models 



For a powerlaw loop length spectrum with a negative spectral index, most loops will 
be small. This is a disadvantage in a lensing search because of image separation: lensing 
events from small loops are likely to be lost due to resolution limits. With CLASS/ JVAS, we 
expect to see few lensing events in the survey over most of the range of plausible powerlaw 
indices and loop density parameters. LOFAR will have the ability to rule out a slightly 
larger portion of the parameter space, for Qioop ^ 10"^ and 7 > —1. When 7 < 0, the 
sky fraction is strongly dependent on the instrument's ability to resolve the smallest loops; 
both CLASS/ JVAS and LOFAR lack the resolution to be effective in this regime. With 
the high resolution of SKA, however, nearly all loops will be resolved for powerlaw indices 
—3 < 7 < with L^, = 1 kpc, and the large number of compact radio sources available 
will allow us to rule out to 3-a a density parameter in loops of ~ 10"^ if no loop lensing 
events are observed. The results for the powerlaw loop spectrum are presented in Figure 
m (Note that these results could easily be translated into constraints on any population of 
point masses with a powerlaw mass spectrum, which may include some models of primordial 



black holes 



TABLE L Properties of Gravitational Lensing Surveys 



Survey 


Resolution 


Number of Sources 




(arcsec) 


(observed/expected) 


CLASS/ JVAS 


0.3 


11685 


Ext. LOFAR 


0.8 


1.5 X 10^ 


SKA 


0.01 


108 
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FIG. 2: Contours for 10 and 100 expected events using the powerlaw loop spectrum model for 
CLASS/JVAS (C), LOFAR (L) and SKA (S), for a source population at redshift z=1.2. The x- 
axis is the powerlaw index 7 (see Section [II B ip and the y-axis is the density parameter in loops, 
^loop- 10 expected events corresponds to a constraint at 3a if no events are observed. 

2. One-scale model. 

Our results for a one-scale model spectrum of loops are presented in Figure [3] for the 
CLASS/JVAS, LOFAR and SKA surveys. The one-scale model spectrum has two limiting 
cases of interest, depending on the relationship between Gfi/c^ and a. 

When a < IGG/i/c^, the upper cut-off in the spectrum is below the peak, and thus the 
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spectrum is approximately a powerlaw with spectral index +2 and is dominated by loops 
near the upper cut-off at 2alHo- -^^^ ^ powerlaw spectrum as in Equation H] with 7 > 0, the 
mass density goes as 

with A^o ~ ^loopl {[GlJ'/ c^]Lmax)- This means that the mass density p ~ ^loop ~ Ija. 
Ignoring resolution effects, the sky fraction is proportional to the mass density, since 9\ ~ 
Mioop, so in this regime, the sky fraction is independent of Gfi/c^. When a is small, the 
resolution of the instrument also decreases the sky fraction because the smallest image 
separations will be lost. 

When a » Gfi/c^, the mass in loops is dominated by the region of the spectrum 
with a spectral index of 0. For 7 = 0, p ~ {G fi/c^)No In (Lmax/Lmin)- Since the low- 
mass end of the spectrum gives a small contribution, one can approximate the lower cut-off 
as the peak position: L^m ~ Lp^ak ~ Gji/c^. Lmax is the upper cut-off which goes as 
a. In this case, Nq ~ Qioop/ {[Gfi/c^]\na/[Gfi/c^]). Since in the limit that a » Gfi/c^, 
^loop ~ {Gfi/c'^)\na/{Gfi/c'^), one finds that p ~ {Gp/c^)\na/ (Gp/c^). In this regime, the 
sky fraction depends only very weakly on a. 

If the loop spectrum follows the one-scale model, we can rule out regions of the parameter 
space to 3-cr even with the CLASS/ JVAS survey, as seen in Figure [31 This plot also shows 
that a large region of parameter space can be ruled out with LOFAR or SKA. With SKA, 
a 3-0" constraint will be achievable down to Gp/c^ ~ 10~^ for large values of a. In the next 
section, we discuss how loop lensing surveys can be more powerful and more direct than 
current methods to constrain cosmic strings with pulsar timing. 



B. Comparison with other constraints 

Currently, the strongest limits on cosmic strings come from constraints on the stochastic 
gravitational wave background obtained from the timing of millisecond pulsars j^J]. In the 
context of the one-scale model, pulsar timing constrains a combination of Gp/c^ and a. 



In Figure [3l we include the pulsar timing constraint from [45| for two different values of 
the density parameter in gravitational waves for comparison with the limits we can obtain 
with radio surveys. The limit obtainable from SKA can improve upon the pulsar timing 
constraint by up to two orders of magnitude in Gp/c^ in some regions of the parameter 
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FIG. 3: Contours for 10 and 100 expected events using the one-scale loop spectrum model for 
CLASS/JVAS (C), LOFAR (L), and SKA (S), for a source population at redshift z=1.2. We 



also plot the pulsar timing constraint from 



451] (black dottec 



r^g/i^ < 2 X 10® and the lower is for ^Igh'^ < 2 x 10^ - see 



lines). The upper dotted line is for 



45l | for a discussion of the uncertainty. 



In this plot, the x-axis is the one-scale model parameter a (see Section [11 B 2\i and the y-axis is the 
dimensionless loop tension Gn/c?. 



space. 

However, a simple comparison between limits from the gravitational wave background 
and from loop lensing surveys can obscure the most important feature of a lensing survey, 
which is that it is a direct search. While constraints from pulsar timing rely on models 
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of the loop radiation and simulations of the behavior of loop networks, lensing relies only 
on the presence of a mass concentration, the effects of which are directly observed in the 
image distributions. Any model of a cosmic string population will make simple and testable 
predictions for gravitational lensing surveys. 



IV. DISCUSSION 



We have found that searching for lensed compact radio sources in future surveys can 
greatly improve constraints on cosmic string networks. Extending the excluded region of 
parameter space, or alternatively observing a cosmic string, will allow us to better under- 
stand physics beyond the Standard Model. We have shown that if no cosmic string loop 
lensing events are observed in upcoming radio surveys, much of the currently interesting 
parameter space can be ruled out, but further work is needed to determine what one could 
learn about the nature of cosmic strings if a lensing event were actually to be observed. For 
the sake of obtaining a constraint, it is sufficient to discuss the scenario in which all lensing 
events in a survey can be accounted for by the detection of a conventional lensing object 
(such as a galaxy or cluster) in follow-up observations, so there is no need to invoke cosmic 
strings. In all known strong gravitational lensing events to date, observations have identified 
the source of the lensing potential as a standard baryonic or dark-matter structure. 

In this work, we have tried to make conservative assumptions in order to estimate the 
lensing signal, but including additional effects is only likely to improve our limits. For in- 



stance, the statistically complete sample of lenses in CLASS [36| includes only those systems 
with fiux ratio less than 10:1, but future radio interferometers with higher dynamic range 
may be expected to easily improve on this ratio by at least an order of magnitude. This 
would significantly increase the effective lensing area for each loop to as much as an order 
of magnitude above the area within the Einstein radius. 

In any survey, brighter events are more likely to be above the fiux limit, and so lensing 
surveys will see a disproportionate number of high-magnification events. By not explicitly 
including the effect of this magnification bias in our estimate, we are underestimating the 
number of lensed sources. 

The detailed properties of the cosmic strings under consideration will also likely enhance 
their prospects for detection. One such property is the reconnection probability P, which 
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is unity for ordinary cosmic strings 46|], but can be significantly less than unity for cosmic 
strings originating from superstring models 
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the string number density increases as ~ p-o.e which would proportionally increase the 
number of expected lensing events. In addition, there may be several distinct string popu- 
lations present: for example, in brane inflation scenarios multiple types of cosmic strings, 
labeled by coprime integers (p, q), would be produced. The total number of strings would be 



49| . enhancing the lensing signal 



enhanced by the number of different populations present 
by the same factor. 

Further work on the statistical properties of loop lensing will help to sharpen our con- 
straint. We expect that considering lensing by loops in more general configurations (includ- 
ing oscillations) will only slightly affect the total lensing cross section. However, understand- 
ing the effects of configuration may be important for determining the detailed properties of 
individual loops detected in lensing surveys. 

Experiments such as LISA will revolutionize the observation of gravitational waves, in 
particular the stochastic background, to which cosmic strings are expected to contribute 
a distinctive signature 50|]. In comparison with limits from millisecond pulsars, LISA is 
expected to push back the minimum detectable string threshold by ~ 7 orders of magnitude 



44( 1 . Although these tensions are inaccessible to the surveys we mention in this work, a 



lensing search is a complementary technique in that in addition to requiring fewer assump- 
tions about the behavior of strings, it involves a local phenomenon rather than a spectral 
signature from a population. 

We hope that our results will motivate the optimization of future surveys for events of 
this nature. The next generation of radio surveys have the potential to directly detect cosmic 
string loops, opening a new window into the physics of the early universe. 
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APPENDIX A: ONE-SCALE MODEL LOOP SPECTRUM 

In this Appendix we describe the process by which the loop length distribution ([6]) and 
total mass in loops ([7]) are derived in the context of the one-scale model (OSM). The equa- 
tions and parameter values we use are those of jj]. According to the OSM, the energy 
density in long strings scales like where I is the particle horizon, defined by 

i{t) = a{t) . / 4^ = 3rf, (Al) 



JO a{t') 

and the last equality holds during matter domination. Applying energy conservation then 
determines the rate at which the long string network sheds energy in the form of loops. 
The loops formed at time t are all taken to have an initial length al{t), where a is a 
dimensionless constant. These assumptions imply that the absolute number of loops 
created in a comoving cubical volume of size i? on a side is given by 

dTV _ a{tfR^ C^_^iHo .. . 

d£ £^ ' a Aa P ^ ^ 

where C = 5.3 during matter domination. To obtain the second equality, we have taken 
a = 1 at a fiducial time to, and used a(t) = {t/toy^^ = (i/2iHo)'^^^, where ino = c/H(t) is 
the Hubble length at to- We have further taken R = ino, so N is the number of loops in a 
comoving (cubical) horizon volume. 

Once loops have formed, they begin to shrink as they lose energy by emitting gravitational 
radiation. The energy loss rate is 

d E 

^ = -rG/i^c, (A3) 

where F ~ 50 is taken from simulations. This means that the length of a loop that was 
"born" at time is 

Lit,tB)=fraiitB)-^{t-tB), (A4) 

where /r ~ 0.7 is a factor accounting for the shrinking of the loop as its initial relativistic 
velocity is redshifted by Hubble expansion. During matter domination the time at which 
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a loop "dies" (completely evaporates) is related to its time of birth by 

tD 



TGfi 



(A5) 



This also means that a loop of length L, observed when the particle horizon is i, was born 
when the particle horizon was a size given by 



3p2 



(A6) 



In our universe, if /? > 5 x 10"^, then all the loops created during radiation domination will 
have evaporated by the present day. We will assume this is so in the following. 

Our first goal is to calculate the spectrum of loop lengths, evaluated at the fiducial time 
to- We proceed by writing 



AN 



dN dis 



(AT) 



to 



diB dL 

where is the particle horizon length when a loop currently of length L was created. The 
first derivative on the right hand side is given by (lA2p . and the second is obtained from 
(]A6p . Combining them, we find 



dA^ 
'dT 



Cc^ 



L/i 



m 



to 
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2 ' 



(A8) 



which, once the values for fr,T and C are substituted, becomes precisely equation (I6j). 

To compute the density parameter in loops, we first calculate the total length in loops. 
This requires the integral identity 



(1 



AL/B 



dL = AB 



L/Bf 

which applies to our case with the choice. 



In 1 



Ljf 

'b 



LJB 
l + L,/B 



(A9) 



(AlO) 



3c2 

The upper cutoff is necessary because the largest possible loops are those being created right 
at to- The integral formula gives 



- (tot) 
''loop 



3/. 



-y 



In ( 1 + -y 



(All) 



where we have taken y = a/^Gfi/c^). From the standard definition of the density parameter, 
we have 



'-HO 



Finally, substituting the values of fr, T and C we obtain equation ([7]). 
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